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EXPONENTIAL CONVERGENCE TO EQUILIBRIUM IN A 
POISSON-NERNST-PLANCK-TYPE SYSTEM WITH NONLINEAR 

DIFFUSION 


Abstract. We investigate a Poisson-Nernst-Planck type system in three spa¬ 
tial dimensions where the strength of the electric drift depends on a possibly 
small parameter and the particles are assumed to diffuse quadratically. On 
grounds of the global existence result proved by Kinderlehrer, Monsaingeon 
and Xu (2015) using the formal Wasserstein gradient flow structure of the sys¬ 
tem, we analyse the long-time behaviour of weak solutions. We prove under 
the assumption of uniform convexity of the external drift potentials that the 
system possesses a unique steady state. If the strength of the electric drift is 
sufficiently small, we show convergence of solutions to the respective steady 
state at an exponential rate using entropy-dissipation methods. 
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1. Introduction 

We are concerned with the long-time behaviour of the following parabolic system 
in three spatial dimensions: 


(1) 

dtu = div(uD(2M + U + eip)), 

(2) 

dtv = div(uD(2u + V — ei/’)), 

( 3 ) 

—A-ip = u — V, 


for nonnegative u,v : [0, oo) x K.^ —>■ [0,oo], together with an initial condition 
m(0, ■) = > 0, t( 0, •) = > 0 on Equations (1)&(2) are coupled by means 

of Poisson’s equation (3) giving 


Ip = G * {u — v), 

with Newton’s potential G in i.e. 

G(a;) = for x ^ 0. 

A'K\x\ 
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For the external confinement potentials U,V G C^(K^), we assume that they grow 
quadratically: 

Boundedness of second derivatives: ||D^17 ||l“ < oo, ||D^{7 ||loo < oo, 
Uniform convexity: > AqI, D^U > AqI, 

in the sense of symmetric matrices, for some Aq > 0. We may assume without loss 
of generality that U,V >0. Finally, e > 0 is a fixed parameter. 

System (l)-(3) possesses a gradient flow structure on the space X = 3^2 x 3^2, 
where i ^2 denotes the space of absolutely continuous probability measures on 
with finite second moment, endowed with the metric 

d((u,M), (u,u)) := ^/W2{u, uY + W2{v, v)^, 

where W 2 is the L^-Wasserstein metric. The corresponding energy functional £ : 
X —> K U {+ 00 } reads 

+ uC/+ uU + ||Di/)p) dx, if (u,'!;) e X 

£{u,v) := , 

1 + 00 , otherwise. 

It has been shown by Kinderlehrer et al. [19] that, given (u°, u°) € Xfl x L^), 
a global-in-time weak solution (u,u) : [0,oo) —>■ X to (l)-(3) exists and can be 
constructed as the continuous-time limit r —>■ 0 (in a sense to be specified below) 
of the minimizing movement scheme with step size r > 0: 

(u?,u?):= (u°,u0), 

(u", u”) e argmin (—d(('u, u), -I-f (u, u)) for n € N. 

(u.«)ex ' 

We will summarize important properties of the sequences (u”,u")„gN and their 
limit (u, v) in Section 2 below. 

In this work, we are interested in the behaviour of the aforementioned weak solu¬ 
tion {u,v) to system (l)-(3) as t —>■ 00 . First, we characterize the set of equilibria: 

Theorem 1.1 (Existence and uniqueness of stationary states). For every e > 0, 
there exists a unique minimizer {uoo,Voo) G x lU^’^) of £ on X. (rtoc'^oo) is 

a stationary solution to (l)-(3) and satisfies 

(fi) U 

(6) Uoo = ]^\Cv - u -f e'0oo] + , 

f’oo ■- Gr ^ {Uqq Uoo); 

where Cu,Cy G M are such that ||uoo||li = 1 = ||uoo||li; [-j-i- denoting the positive 
part. For every a G (0,1), Uoo,Voc G (7°’“ with compact support and ip G 

Second, we prove for sufficiently small coupling strength e > 0 exponential con¬ 
vergence to (Uoo,Voo)- 

Theorem 1.2 (Exponential convergence to equilibrium). There are constants e > 0 
and L > 0 such that for all A > 0, there exists Cs > 0 such that the following is 
true for every e G (0,e) and arbitrary initial conditions (u^,v^) G Xfl (L^ x L^): 
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The weak solution (u^v) to (l)-(3) obtained as a limit of the scheme (4) converges 
to (uoo,Voo) exponentially fast with rate := Ao — Ae > 0 in the following sense: 

W 2 {u{t,-),Uoo) +'W 2 {v{t,-),Voo) + •) “ Uoo||l 2 + ||^^(^, •) “ Voo\\l^ 

(7) 

< Cs (£(uo, Vo) - f (woo, foo) + for all t > 0. 

System (l)-(3) may arise as a model for the dynamics of a system consisting of 
positively and negatively charged particles (e.g. ions) inside some electrically neu¬ 
tral surrounding medium (e.g. air, water). For further details on the mathematical 
modelling of those phenomena, we refer to the monographs [23, 18]. Here, both 
species are confined by means of external potentials U and V and are assumed to 
diffuse nonlinearily - with a diffusive mobility depending linearily on the concentra¬ 
tions u and V, respectively. We assume the Poisson coupling by means of equation 
(3) to be suitably weak (e <C 1), i.e. the drift induced by electromagnetic force to 
be small. The quantity 1 corresponds to a large relative permittivity (dielec¬ 

tric constant) of the surrounding medium. A similar system has been considered 
by Biler, Dolbeault and Markowich [7]. There, a time-dependent coupling e(t) was 
introduced, with the crucial assumption that e(<) —>■ 0 as t —>■ oo, i.e. asymptot¬ 
ical damping of the electrostatic potential. Under relatively general requirements 
on spatial dimension, external potential and diffusive nonlinearity, convergence to 
equilibrium as t ^ oo is proved for sufficiently regular solutions. Here, we do not 
require asymptotical damping of the Poisson coupling, that is, the system at hand 
still constitutes a coupled system even in the large-time limit t —>• oo. To the best 
of our knowledge, our rigorous result on exponential convergence of weak solutions 
is novel in the case of genuinely nonlinear diffusion on multiple space dimensions, 
even in the small coupling regime e -C 1. Partial results have been obtained in one 
spatial dimension [15] or for space-dependent diffusion [5] only. 

In contrast to that, the case of linear diffusion has already been treated almost 
exhaustively. In the articles [3, 6, 4] preceding [7], it was shown that the rate of ex¬ 
ponential convergence to equilibrium of the system without coupling, for uniformly 
convex potentials, is (almost) retained for coupled systems. There, the strategy of 
proof is mainly based on applications of generalized Sobolev inequalities the deriva¬ 
tion of which require the use of a Holley-Stroock-type perturbation lemma. Seem¬ 
ingly, such a strategy might not be applicable in the setting of nonlinear diffusion. 
On the other hand, systems of the form above possess (at least formally) a gradient 
flow structure (w.r.t. e.g. the T^-Wasserstein distance) which also is of use for the 
analysis of the system - and, in contrast, does not at all require linear diffusion. 

Variational techniques related to gradient flows in (transportation) metric spaces 
[27, 2] have been applied to a variety of evolution equations, e.g. in [17, II, 26, 13, 
14, 12, I, 16, 24]. The variational minimizing movement scheme [17] provides a 
key tool, in combination with generalized convexity assumptions on the respective 
free energy or driving entropy functionals [25], for the investigation of existence 
and long-time behaviour of solutions to nonlinear evolution equations with gradient 
structure. Recently, also genuine systems of equations were object of study in this 
context. Using the minimizing movement scheme on an appropriate metric space, 
existence of weak solutions has been proved in several cases, e.g. for Keller-Segel- 
type systems [8, 9, 10, 28, 30, 29] or others [20, 21, 31, 19]. Using the minimizing 
movement scheme to obtain convergence to equilibrium is rather novel in the case of 
genuine systems. The method applied here has first been used for Keller-Segel-type 
models in [30, 29] leading to similar results as Theorem 1.2 here. 
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The basis of our strategy is the fact that the uncoupled system {e = 0) defines a 
Ao-contractive flow, since £ then is Ao-convex along geodesics in (X, d). However, 
geodesic convexity of £ is lost if e > 0. Still, one can prove (see [19]) that a 
continuous flow in X of system (l)-(3) exists and admits a unique steady state 
{uoo,Voo) (see our Theorem 1.1). The cornerstone of the proof of Theorem 1.2 is 
the existence of an auxiliary functional £ which is “close” to £ (for both e > 0 
and e = 0): First, £ is Ag-convex along geodesics in X, with a slightly smaller 
convexity modulus 0 < A^ < Aq, and is decoupled in its arguments u and v. Hence, 
known results on gradient flows for scalar porous-medium type equations apply for 
the auxiliary gradient flow of £. This auxiliary flow can be used - with the 
almost classical flow interchange technique from [24] - to estimate the dissipation 
of £ along the continuous flow given by the free energy £. We seek to eventually 
apply Gronwall’s lemma for £. Since £ and £ differ by a “small” - but non-convex 
- functional, cross-terms occur in the entropy-dissipation estimate and have to be 
controlled by suitable a priori estimates. For small coupling strength, we arrive in 
the end at an exponential estimate with an again smaller rate 0 < < Ag < Aq. 

Clearly, this strategy requires Aq > 0, i.e. uniform convexity of the external po¬ 
tentials U and V, which is not needed for proving existence (see [19]). As in [30], we 
deal with quadratic diffusion only since the right dissipation estimates do not seem 
to be at hand in the general case. One last comment is due about the scaling of our 
exponential estimate (7) in Theorem 1.2: For uniformly contractive gradient flows, 
one expects the difference of initial and final energy to enter the estimate with a 
power 1/2, corresponding to 6 = —1/2. However, due to non-convexity of the free 
energy £, only positive S can be obtained with our strategy here. Nevertheless, the 
initial condition only appears via its energy. 

This paper is organized as follows: First, we recall general facts and definitions 
for gradient flows in metric spaces and a result on the global existence of solu¬ 
tions to the system at hand. In Section 3, we prove Theorem 1.1 on existence and 
uniqueness of steady states. Section 4 is devoted to the introduction and investi¬ 
gation of the auxiliary entropy mentioned above. There, we also derive a central 
entropy-dissipation estimate for our forthcoming analysis, using the flow interchange 
technique. Finally, Theorem 1.2 is proved in Section 5. There, we first derive an 
additional a priori estimate on the auxiliary entropy holding for large times. In 
consequence, exponential convergence is proved. 


2. Preliminaries 


2.1. Geodesic convexity and gradient flo-ws. In this section, we will briefly 
mention relevant definitions and facts on gradient flows in metric spaces. For a 
more thorough presentation, we refer to [2, 27]. 

Throughout this paper, D and denote the spatial gradient and Hessian, re¬ 
spectively. By abuse of notation, we often identify an absolutely continuous measure 
with its Lebesgue density. A sequence {pn)n&n of probability measures on is said 
to converge narrowly to some limit probability measure fi if for all continuous and 
bounded maps />: —>■ R, one has 



cj){x) dp,{x). 
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For the metric space (^ 2 ,W 2 ), the following is true: A sequence in 

converges w.r.t. the L^-Wasserstein distance W 2 if and only if both fin ^ nar¬ 
rowly and the sequence of second moments converges: 

\x'^ dp{x) for p S 3^2- 

A functional ^ : X —>■ M U { 00 } defined on some metric space {X, d) is called 
X-geodesically convex for some A G R if for every wo,Wi € X and s € [0,1], one has 

-A(ws) < (1 - s)Al(wo) + sAl(rci) - -^5(1 - s)d^{wo, wi), 

where Wg : [0,1] —>• X, s 1 —>■ uig is a geodesic curve connecting wq and wi. We recall 
two important classes of A-convex functionals (see e.g. [2, Ch. 9.3], [27, Thm. 5.15]) 
on the space (.^ 2 , W 2 ): 

Theorem 2.1 (Criteria for geodesic convexity on {3^2, XX 2 )). The following state¬ 
ments are true: 

(a) Let h G ^^([O, 00 )) be given, and define a functional A on 3^2 by Al(w) := 
/jjd h{w{x)) dx for w G £^2 D L^. If h{0) = 0 and r 1 —>■ r'^h{r~‘^) is convex and 
nonincreasing on (0,oo), A is 0-geodesically convex and lower semicontinuous 
in ('^ 2 ,W 2 ). 

(b) Let a function W G C°(R‘^) be given, and define a functional A{p) := W dp 
for all p G 3^2- If W is X-convex for some A G R., .A is X-geodesically convex 
in {3^2, W 2 ). 

As to the notion of gradient flow, we use the following characterization: 

Definition 2.2 («:-contractive flow). Let Al : W —>■ RU { 00 } be a lower semicontin¬ 
uous functional on the metric space {X, d). A continuous semigroup on {X, d) 
is called n-flow for some k G M, if the evolution variational estimate 

+ '^d^i^fM,w)-£AiSf{w)) < A{w) 

holds for arbitrary w, w in the domain of A, and for all t > 0. 

We recall some facts on gradient flows of convex functionals on £^2 ■ 

Theorem 2.3 (Gradient flows of geodesically convex functionals on {£^ 2 , W 2 ) [2]). 
Let A : £1^2 —>■ R U { 00 } be lower semicontinuous and X-geodesically convex w.r.t. 
the distance W 2 . The following statements hold: 

(a) There exists a unique n-flow, with k := X, for A. Its corresponding evolution 
equation can be written as 

dtSfiw) = div (^Sj^(u;)D , 

if A is sufficiently regular. There, stands for the usual first variation of the 
functional A on L^. 

(b) There exists exactly one minimizer Wmin of A, for which the following holds: 

A„, 2 , , ... ., , 1 A(w) — A(S-^(w)) 

(8) TtWKw, Wmin) < A{w) - Al(lCmin) < 1™ -T-• 

Z Za h\0 fl 

One of the cornerstones of our analysis below is the following theorem: 


lim m. 2 {pn) = Tna. 2 {p), with m. 2 {p) := / 

n-i>oo 
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Theorem 2.4 (Flow interchange lemma [24, Thin. 3.2]). Let B be a proper, lower 
semicontinuous and X-geodesically convex functional on {X,d). Let furthermore A 
be another proper, lower semicontinuous functional on {X, d) such that Dom(y4) C 
Dom(S). Assume that, for arbitrary r > 0 and w G X, the functional ^d{-,w)'^ +A 
possesses a minimizer w on X. Then, the following holds: 

B{w) + tT)^A{w) + ^d^{w, w) < B(w). 

There, Id’^A{w) denotes the dissipation of the functional A along the gradient flow 
S® of the functional B, i.e. 


T>’^A{w) := limsup 


-4(w) -.4(Sf(w)) 
h 


2.2. Minimizing movement and existence of solutions. In this subsection, 
we recall the results proved by Kinderlehrer et al. in [19] in our specific setting. 

Proposition 1 (Minimizing movement [19, Prop. 3.3]). Let r > 0 and (u°,n°) € 
Xn (L^ X L^) be given. Then, the sequence (w",i;")„gN defined by the minimizing 
movement scheme (4) is well-defined with (u",i;") G Xfl x for all 

n G N. By definition, the sequence (£’(«", i;"))„gN is nonincreasing. 


Define for r > 0 the discrete solution (ut, Vr) '■ [0, oo) —>■ X by piecewise constant 
interpolation, that is 

{Ur,Vr){,0) := {u°,V°), 

{Ur, r'r)(0 •= (^r I ) fo'' ^ ^ ((^ ~ I)''") n > 1. 

The following main result of [19] about the existence of solutions to (l)-(3) is at 
the basis of our subsequent analysis: 


Theorem 2.5 (Existence of solutions [19, Thin. 2]). Let £ > 0 and U,V as men¬ 
tioned above be given. Define, for initial conditions {u^, v^) G Xn(L^ x L^) and each 
r > 0 a discrete solution {ut,Vt) by (4)0(9). Then, there exists a sequence Tk\0 
and a map {u,v) : [0,oo) x —>• [0,oo]^ such that for each t > 0, Ur,^{t) u(t) 
and Vri,{t) v{t), both narrowly in as k ^ oo. Moreover, {u,v) is a solution 
to (l)-(3) in the sense of distributions, it attains the initial condition and one has 
for each T>0: 

u,v G C^^^{[0,T]-{3^2,W2)) n L^{[0,T]-L'^) n L\[0,T];W^^^), 
£{u{T),v{T))<£{u^y). 


3. The equilibrium state 
In this section, we prove Theorem 1.1. 

Proof. Existence: Trivially, £ is bounded from below. Hence, there exists a mini¬ 
mizing sequence {uk,Vk)k&i in Xfl (L^ x L^) with lim £{uk,Vk) = inf £{u,v). 

k—>-oo {u,v)£'K 

Thus, we have for some C > 0 that ||Mfe||L 2 < C, llrifejlLS < C for all k G N. 

Moreover, using the Ao-convexity of U and V on one obtains supm 2 (ufc) < oo 

fceN 

and supm 2 (ffe) < oo with the help of the elementary estimates U{x) — > 

fceN 

and V{x) -V{x^iJ > ^\x\^ - (with the unique min- 

imizers of U and V on R^, respectively). We infer with the Prokhorov 
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and Banach-Alaoglu theorems that there exists a subsequence (non-relabelled) and 
a limit (uoo, Vca) G X fl x L^) such that Uk Uoo and Vk Voo both narrowly 
as probability measures and weakly in as /c —>■ oo. With respect to these con¬ 
vergences, £ is lower semicontinuous. In fact, this is obvious for the quadratic and 
linear terms in £ since (7 and V grow quadratically. For the last term containing 
the Dirichlet energy i||D'0|||2, we refer to [19, Prop. 6.1] for a result on lower semi¬ 
continuity w.r.t. weak convergence. Hence, it follows that (uoo,Uoo) is indeed a 
minimizer of £ on X and hence also a steady state of (l)-(3). 

Uniqueness: We claim that £ is uniformly convex with respect to the flat distance 
induced by the product norm jj • which implies the uniqueness of minimizers. 

Indeed, for all (u, v), (u', v') € Xn {L^ x L^) and all s € [0,1], we have, thanks to 


(10) / |D(G + w)pdx= / (G + w)ri;dx= / / w{x)G{x — y)w(]j) dxdy 

JK? Jv? Jv? Jk? 


/R 3 Jv? 

which holds for all u; € X fl that 

ds2 


£(u + s{u — u), V + s{v' — v)) 


s=0 


= ( [2(u' — -I- 2{v' — vY + £{{u' — u) — {v' — v))G * {{u' — u) — {v' — u))] da: 

aR3 


> 2||(u'-M,u'-u)||^2^i2, 


so £ is 2-convex w.r.t. the distance induced by || • Wl'^xL'^- 

Euler-Lagrange equations: Since (uoc'Coo) is the minimizer of f, the following vari¬ 
ational inequality holds: 

£{Uoo + SU,Voc, + sv) 

s=0 

(11) = / [2Uoo + U + sG * {Uoo - Voo)]udx 

JR3 

+ / [2Voc-V - £G*{Uoo-Vao)]vdx, 

JK? 

for all £t, u such that both Uoo -|- ft > 0 and Voo -I- u > 0 on and udx = 
0 = /jg 3 vdx. In order to prove (5), we set v := 0. Let (() : —>• M be such that 

/r3 (pdx <\ and (j) -I- Uoo > 0 on R^. The choice 

4>dx 

is admissible for u in (11), hence (recall our notation := G * (uoo — Uoo)) 




( 12 ) 


0 < 


I (2Uoo d~ U £'^oc 
JR3 


Cu)4>dx, 


with 

Cu := / (2u^ + Uu ) dec € R. 

Jr3 

If Moo (a:) > 0 for some x G R^, we are able to choose (f) supported on a small 
neighborhood of x and to replace by —(p in (12) and obtain 

2u^{x) + U{x) + eipocix) = Cu- 
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If Uooix) = 0 for some x, one has U{x) — eipooix) — Cu >0, and hence (5) is true in 
both cases. The equation for Voo (6) can be derived in analogy. 

Properties: First, since (Moo,foo) are admissible as starting condition (for 

arbitrary r > 0) in scheme (4), we obtain thanks to the minimizing property and 
Proposition 1 that (moc'^oo) S x We now show that '0oo G L°°. To this 

end, let a; G and observe at first that 



\uoo{y) - Voo{y)\ 

\x - y \ 


dy < ||uoo 


^oo IIl^ 


/Bi (x) 


'- y ? 



= 2v^||Moo - J'ooIIls, 


1/2 


independent of x, by Holder’s inequality and the transformation theorem. Further¬ 
more, since \x — y\>l\i y ^ Bi(a:) and llitaDllLi = 1 = IIuooHli, we get 
f \uoo{y) - Voo{y)\ 

iR3\Bi(a:) 


' - y \ 


■ dy < ||uoo -'I’ooIIli sup |a;-y| ^<2. 


Putting both parts together, we see that sup |'0oo(a:’)| < oo. In view of (5)&(6), 

ipoo G L°° implies that Uoo and Voo have compact support since U and V grow 
quadratically as |a;| —>■ oo. By classical results on solutions to Poisson’s equation 
[22, Thm. 10.2], we then infer that tpoo G (7°’“ for all a G (0,1), since by the 
Gagliardo-Nirenberg-Sobolev inequality, one has {uoc,Voc>) G x L^. Hence, using 
(5)&(6) again, we conclude that Uoo and Voo also are Holder continuous. By elliptic 
regularity for Poisson’s kernel [22, Thm. 10.3], it follows that tp^ € □ 


4. Auxiliary entropy and dissipation 

In this section, we define a suitable geodesically convex auxiliary entropy £ and 
derive the dissipation of the driving entropy £ along the gradient flow S'^ of £. 


Let £ : X —> M U {oo} be defined via 

7r3 + + uoo)u + (v- voo)v 

+e(u - Uoo)lpoc - £(v - Voo)lpoo ] dx, 
if (u, v) € X L^, 

^-boo, otherwise. 

Obviously, £ is proper and lower semicontinuous on (X,d). 


£(u, v) := 


Proposition 2 (Properties of £). There exists Eq > 0 such that for all e G (0,Eo): 
the following statements hold: 

(a) There exists L > 0 such that £ is Xe-geodesically convex w.r.t. d, where Xe := 
Aq — Ts ^ 0. 

(b) The following holds for all {u,v) G Xfl (bF^’^ x bF^’^).- 

||m- Moolis -b \\v-Voo\\l2 

(^^ 3 ^ <£iu,v) 

< [ [u|D(2u-b 1/-b eV’oo)P + i^|E)(2u-b F - eV’oo)!^] dx. 

(c) There exists a constant K > 0 independent of e such that for all (u,v) G X.' 

(14) £{u,v) < £{u,v) - £{uac,Vao) + Ke. 
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Proof, (a) In view of Theorem 2.1, as C is decoupled in its arguments u and v, 
it suffices to prove that there exists C > 0 such that ||D^'!/'oo||l“ < C for all 
sufficiently small e > 0. Let R > 0 such that suppwoo U supp?;oo C 
Since ipoc G thanks to Theorem 1.1, we have sup \dxidxj'f’oo{x)\ < oo 

xGMr+i{0) 

for each pair {i,j) € {1,2,3}^. Consider now x ^ ]B/j+i(0). One easily obtains 
for z 0 that 


d,,d,^.G(z) = 


47r|i 


SziZj 


where Sij denotes Kronecker’s delta. So, using a linear transformation, 


\dxi9xj'4^0o{_^'}\ - 


/ dz,dz.G(z)(uoo(x - z) - Voc{x- z))dz 
JMr(x) 


< - V. 


oo ”oo \\L°° 5 


since for all 2 € Mji{x), one has \z\ > 1 by definition of x. Hence, the desired 
uniform estimate on D'^'ipoo is proved. 

(b) The upper estimate is a straightforward consequence of A^-convexity of C and 
the structure of its Wasserstein subdifferential w.r.t. u and v, respectively (see 
e.g. [2, Lemma 10.4.1]), in combination with (8). For the lower estimate, we 
observe that 


£(m, v) 

' I [ ^oo) “t“ ip "Coo) F {u 'Uoo)(2Woo U -\- 
+ (w - 'Coo)(2z;oo + V - eipoo) ] dx. 


We prove that J^ 3 {u — Uoo)i‘2uoo+ U + eipoo) dx > 0. Since the last term above 
can be treated in the same way, the claim then follows. Using (5), we obtain 


/ (it —'«oo)(2'«oo + U + e-i/ioo) dx 

/ (^ ^oc'jGu dx F / 

J {Cu-U-ei)^>0} 

— Gu I (it itoo) dx F / 
dR3 JI 


u(U F eifoo) dx 


'{C„-f7-£j/'=o<0} 


it(C/ F ef; OO Cu) dx 


> 0 , 


since it and itoo have equal mass (hence the first term is equal to zero) and the 
integrand of the second integral is nonnegative on the domain of integration, 
(c) One has for all (it, x) G X 0 (L^ x L^): 


-{£{U,V) - £{U,V) + £{Uoo,Voo)) 
£ 


/R3 


(it ltoo)V^OO (x ltoo)V^OO 2 !^'^! 


dx 


f 11 

< / ■*/'oo(u-^'- xMoo F -Xoo)dx < 3 ||V'oo||l° 

jR3 ^ ^ 

<K. 
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thanks to (10) and Theorem 1.1. 


□ 


(15) 


According to Theorem 2.3(a), the Ag-contractive flow S'^ =: iU^V) is character¬ 
ized by 

dJA = div [WD(2W + U + e ^^)], 
dsV = div [VD(2V -b - et/joo)] ■ 

Now, we derive the central a priori estimate on the discrete solution: 

Proposition 3 (Dissipation of E along S'^). Let r > 0 and let (w",i;")„gN be the 
sequence defined via the minimizing movement scheme (4). Then, for all n € N.' 

(16) £«,<) + r2?«,<) < <-i), 


(17) 


the dissipation being given by 

T>{u,v) := (l - ^ [ (u\D{2u +U + + v\D{2v+ V - dx 

^ 2 / 

/ {u + v)\I){'ip-'ip^)\‘^dx. 

2 JR3 

Proof. To justify the calculations below, we regularize the flow given by (15). De¬ 
fine, for > 0 and {u,v) € X fl {L^ x L^) the regularized functional 

£i,(u, v) := C{u, v) + vLtiu) + v'H{v), 

with Boltzmann’s entropy 'H{w) := /j^a wlogwdx, which is finite on (cf. e.g. 

[28, Lemma 5.3]). Furthermore, by Theorem 2.3(a), H is 0-geodesically convex on 
^ 2 , so Li, is Ae-geodesically convex w.r.t. d and the associated evolution equation 
to its Ae-flow {U, V) is the strictly parabolic, decoupled system 

dsU = izAU + div [UB{2U + U + efi ^)], 
d,V = JzAV-bdiv[VD(2V-bP-ei/^oo)]. 

Let {u,v) G Xn (LF^’^ x VF^’^). At least for small s > 0, system (18) has a smooth 
and nonnegative solution {U,V) such that (Z/f(s), V(s)) —>■ {u,v) both strongly in 
X and d, as well as weakly in LF^’^ x LF^’^, for s \ 0. Moreover, this local 
flow can be identified with the Ae-flow associated to (see e.g. [2, Thm. 11.2.8]). 
Then, writing 'L := G * ilA — V) for brevity: 

—(W, V) = - / [2UPU + £T]div [vDU + UP)[2U + U + da: 

ds Jr3 

- [ [2V + V- e«']div [izDV + VD(2V -b P - e^/>oo)] da;. 

dR3 

We first focus on the viscosity terms and obtain, using that {U,V) G X: 

- I {[2U + U + sd>]AU -b [2V -b P - £4-]AV) da; 

dR3 

= / {2\DU\^+ 2\DV\^-UAU-VAV - e{U-V)A'i>)dx 

JR3 

= 2\\DU\\l2 + 2||DV||i2 - j {UAU -b VAV) da; -b e\\U - V\\l^ 
>-\\AUU^-\\AVU^. 
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The remaining terms can be rewritten as 


- [ [2U + U + e^-ldiv [KD{2U + U + dx 
- I [2V + V- e^-Jcliv [VD(2V + F - ei/^oo)] dx 

= / U\D{2U + U + dx + [ V|D(2V + 1/- ei/^oo)^ da; 

aR3 Jr3 

+ £ / WD(2W + C/ + ei/;oo)-D(T-V'oo)da; 

JR3 

- e [ VD(2V + F - ■ D(^ - dx 

JR3 

/ (W|D(2W + C/ + £V'oo)P + V|D(2V + F-£V'oo)P)dx 

^ 3!/ 

^ [ (W + V)|D(iE--V;oo)Pdx, 

^ jR3 


using Young’s inequality in the final step. All in all, we arrive at 
—^^£{U,V) > V{Uy) - p(||AC/|Uoo + ||AY|U=o). 

Observing that the terms appearing in T) are lower semicontinuous w.r.t. the con¬ 
vergence of (W,V) —>• {u,v) above, we obtain after passage to the limits s \ 0 
and p \ 0 that D'^£’(u,t!) > Viu^v). Application of the flow interchange lemma 
(Theorem 2.4) completes the proof of (16). □ 


The remaining task will be to establish appropriate bounds on the dissipation 
T’('«",'u”) in terms of £(■«”, x”) in order to apply a discrete Gronwall lemma and 
to conclude exponential convergence. Note that, in view of (8), it will be enough 
to control the second part of ^(m", w"). 


5. Convergence to equilibrium 

In this section, we complete the proof of Theorem 1.2. Our strategy is as follows: 
First, we derive a uniform bound (independent of £ and the initial condition) on 
the auxiliary entropy C for sufficiently large times. This brings us into position to 
prove a refined estimate on the dissipation T) strong enough to infer exponential 
convergence of C to zero. 

In the following, for r > 0, we denote by (it", ?;")„£« a sequence given by the 
minimizing movement scheme (4). 

5.1. Boundedness of auxiliary entropy. We first need an additional estimate 
for the dissipation terms in (17). 

Lemma 5.1. There exists a constant 6 > 0 such that for all e £ (0,£o) and all 
(ii,i;) G Xn (Wi’2 X 

||m|Ii,3 < (1-I- [ u|D(2it-I-[/-I-£i/)oo)P dx 

V JR 3 

lklll,3 < A-I- [ v\D{2v+ V - dx 

V JR3 


( 19 ) 
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with the convention that the respective right-hand side is equal to +oo if 
u\D{2u + [/ + or z;|D(2i; -\-V — £'i/'oo)P is not integrable. 

Proof. We shall prove the statement for u; the other one can be shown analogously. 
We assume that the r.h.s. is finite. Expanding the square and integrating by parts, 
one has 

/ u|D(2u + [/+ ei/)oo)P da; 

Jk3 

= J + e'0oo) + u|D(C^+ eV'oo)P^ da;. 

Since AU and Aipoo = Voo — Uoo are essentially bounded, we obtain 
^||Dm ^/^|||2 < [ u|D(2u+[/+ e'i/;oo)P dx + C||u||| 2 , 

a JR3 

for some constant C > 0. By the triangle and Young inequalities, one has ||w |||2 < 
2||Wcso|| 1^2 + 2||r/, Woo II 2^2- For small e > 0, we can use (13) and arrive at 

^||Du^/^|||2 < [ fl +u|D(2u +C/+ £'0oo)pdx + 2C'||-«oo||i2. 
y Jr3 \ wy 

On the other hand, with the Lr’-interpolation and Gagliardo-Nirenberg-Sobolev 
inequalities, we have (recall ||u||ii = 1): 

llulUs < \\u\\%%\\]!,^ = ||n3/2||V2 < c"||Du3/2|| 

Raising to the fourth power, we end up with (19). □ 


We now derive a uniform bound on C for large times. 


Proposition 4 (Boundedness of C). (a) There exist ei G (0,£o), E' > 0 and M > 
0 such that for all e G (0, £i), all t > 0 and all n G N.' 

(20) (1 + 2A;t)/:«,O < + rsM, 


where A(. := Aq — L'e > 0. 

(b) Define, with the quantities from (a) and fixed, but arbitrary 5 > 0; 
Mei 


M' := 


Tn := max 0 


2(Ao-L'£i) 

1 + 2(5 


> 0 


and 


2A' 


log 


S{U°,V°) - £{Uao,Voo) + K£1 

W' 


> 0 , 


where K > 0 is the constant from (14). Then, there exists r > 0 such that for 
all £ G (0,£i), r G (0,f] and n G N with nr > To, one has 

(21) £«,<)< 2M'. 


Proof, (a) We first estimate the last term appearing in 'D{u,v) from (17). Let 
{u,v) G X n (W^’^ X W^’^). By Holder’s inequality, 

(22) [ (M + u)|Di/;pdx < (||u||i3/2 + ||i;||^3/2)||Di/;|||6. 
dR3 

The term involving the gradient of i/; can be treated with the Hardy-Littlewood- 
Sobolev inequality (see for example [22, Thin. 4.3] or [19, Lemma 3.1]) which 
is applicable for Poisson’s kernel G: 

(23) llDiAllie < Cl|u - v\\l. < 2C\\u\\l. + 2C\\v\\l., 
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for some constant C > 0. Combining (22)&(23), using HuHli = 1 = ||w||z,i 
again, the L^’-interpolation inequality yields for some /3,/3' £ (0,1): 


/ (m + p)|D('0 —'0oo)p dx < 2 / {u + v)\Dip\‘^ dx + 2 / (u + v)\Dipac\‘^ dx 

Jr3 Jr 3 yR3 

< 4C'(||u||^3||w|li^ + l|M||i3|k|li^ + IMlL3||M||i^ +\\v\\l3\\v\\%) 


+ 2 f {u + n)|D^ooP dx 
JR^ 

<c'(Nli3 + IHIi3 + i), 


for some C > 0, by Young’s inequality and thanks to finiteness of ||Di/!oo||l“- 
Now, we apply (19) and obtain 

V{u,v)> (l-^il + C")) f (n|D(2u + 17 + £V’oo)|" + i^|D(2i; + F-£V’oo)P)dx 
-£M, 


for suitable C" > 0 and M > 0. For £ < we further conclude by (13) 

that 

2?(m, v ) > 2Ae (^1 ~ 2 ~ 


Insertion into (16) yields (a). 

(b) We hrst prove the following explicit estimate for all r > 0 and n £ N U {0} by 
induction over n: 


£«,<)<(£:( 


w°) — £{i 


OO 5 '^OO 

/ ^\—n 


(24) 


Me 


) + Kei){l + 2A(.r) 


+ ^(l-(l + 2A'r)-"). 


Indeed, the claim holds for n = 0 thanks to (14). If it holds for an arbitrary 
n £ N U {0}, we obtain with (20): 

/:«+\<+i) < (l + 2A(,T)-i£«,0 + (l + 2A(,r)-V£M 

< (1 + 2A(.r)“(”+^^ {£{u°, v°) - ^(moo, Poo) + Kei) 

+ ^(1 + 2A;r)-i(l - (1 + 2A;r)-'^) + (1 + 2Kr)-^reM 

= (f (u°, x°) - £-(Poo,Poo) + Ke,)il + 2A't)-('^+i) + ^(1 - (1 + 2A't)-("+i)). 
Let now r > 0 and n £ N with nr >Tq. Thanks to (24), for each (5 > 0, 
£«,0 < ^ + {£{u°,v°) - £{uoo,Voo) + i^£i)exp (^-^log(l + 2 A(,t)^ 

< {£{u°,v°) - £(moo,Poo) + Kei) exp log(l + 2AgT) j + M'. 


Obviously, we obtain (21) in the case £{u°,v°) — £(poo,Poo) + Ksi < M'. 
Consider the converse case. Since lim = 1 there exists s > 0 such 

s-j-O ® 


that 


log(l+s) 


> _±_ 

— 1+25 


for all s £ (0,s]. Henceforth, defining r := ^ yields 
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— TT^ ^ ^ we arrive at the desired estimate by 

definition of Tq: 


£«,<) 

< M' + {£{u°, v°) - £{uoo, Voo) + Ksi) exp 
= 2M'. 


log 


£{u°, v°) - £{uoo,Voo) + 

W' ) 


□ 


5.2. Exponential convergence to equilibrium. We are now able to prove - for 
sufficiently large times - a refined version of Proposition 4(a): 

Proposition 5 (Exponential estimate for C). There exist constants e G (0,£i) and 
L > 0 such that for arbitrary i5 > 0, there exists f > 0 such that for all e G (0,e), 
T G (0, f] and n G N with nr > Tq, we have 

(25) (1 + 2A,t)/:«,<) < £(<-1,<-i), 

with Ae := Ao — > 0 and Tq as in Proposition 4(b)- 


Proof. We write {u,v) instead of for the sake of clarity and consider the 

last term in P(m, v) once more. Using as in the proof of Proposition 4(a) the Holder, 
Hardy-Littlewood-Sobolev and LP-interpolation inequalities (cf. (22)&(23)), we get 
for some C, C" > 0 and /3 G (0,1): 

[ (u + u)|D(^/> -'0oo)P dx 

dR3 

— / ((^ ’^oo) “b (u Uoo) “t“ {Uoo “b ^oo))|U(V^ V^oo)| dx 

< CIKlt - Uoo) - (u - Uoo)||i2 

• {\\u - Uoo||f2||'U- UoolliT^ + ||u - ■!;oo||^2||u - Uoo||^7^ + ||Moo + Voo\\l^/2) 

< C ■ 2C{u, v) ■ C'(l + £(u, v)) < 2CC'{1 + 2M')C{u, v), 


with Young’s inequality, (13) and (21). Now, (25) follows thanks to (16), for suffi¬ 
ciently small £ > 0. □ 


Finally, we prove Theorem 1.2. 

Proof. Consider a vanishing sequence {Tk)kGn such that the corresponding sequence 
of discrete solutions {uT^,,VT,^)keN converges to a weak solution to (l)-(3), in the 
sense of Theorem 2.5. Lower semicontinuity yields for all t > 0: 

C{u{t),v{t)) < liminf £(uT-j,(t),Urs,(t)). By (14) and the monotonicity of £ from 

k—^oo 

Proposition 1, one obtains after passage to fc —>■ oo that 

(26) C{u{t),v{t)) < £{u^,v°) — £{uac,Voo) + Kei Vt > 0. 

Iterating the estimate (25), assuming without loss of generality that G N is suffi¬ 
ciently large, we get in the limit fc —>■ oo that 

(27) C{u(t),v(t)) <2M'exp{-2A^{t-To)) Vt > Tq. 

Actually, (26)&(27) imply that C(u{t),v{t)) < Aexp(—2Aet) for all t > 0, with 
some constant A > 0, the particular structure of which remaining to be identified. 
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Consider the case — £{uoo,Voo) + Kei < M'. Then Tq = 0, so (27) holds 

for all t > 0. In the other case, combining (26)&(27) yields for all t > 0: 

£{u{t), v{t)) < max(£’(M°, v°) — £{uoo,Voo) + Kei,2M') exp(2Aero) exp(—2A£t). 


We insert the definition of Tg and use that A^ < Ag to find 

C{u{t),v{t)) < max((£’(uo,'yo) - Voo) + Kei), 2M') 


f £{uo,Vo) - £iUoc„Vac) + 

■ [ - W' - ) 


Combining both cases yields 

C{u{t), v{t)) < max(£’(Mo, ^ (moo, 'Coo) + ATei, 2M') 


• max 1 


£[uo,vo) -£{u 

OO ; ) + ATei 

W' 


1+2(5 


exp(—2Agt). 


Thus, we can find > 0 such that for alH > 0, the following holds: 

C{u{t), v(t)) < Cs{£{u°, 11 °) - £{uoo,Voo) + exp(-2Aet). 


From here, the desired exponential estimate (7) follows by means of (14) and (8). □ 
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